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Abstract 

We study the twist-3 distribution amplitudes for scalar mesons made up of two valence quarks 
based on QCD sum rules. By choosing the proper correlation functions, we derive the moments of 
the scalar mesons up to the first two order. Making use of these moments, we then calculate the 
first two Gegenbauer coefficients for twist-3 distribution amplitudes of scalar mesons. It is found 
that the second Gegenbauer coefficients of scalar density twist-3 distribution amplitudes for Kq 
and /o mesons are quite close to that for ao, which indicates that the SU(3) symmetry breaking 
effect is tiny here. However, this effect could not be neglected for the forth Gegenbauer coefficients 
of scalar twist-3 distribution amplitudes between oo and /o- Besides, we also observe that the first 
two Gegenbauer coefficients corresponding to the tensor current twist-3 distribution amplitudes for 
all the ao, Kq and /o are very small. The renormalization group evolution of condensates, quark 
masses, decay constants and moments are considered in our calculations. As a byproduct, it is 
found that the masses for isospin 1=1, \ scalar mesons are around 1.27 ~ 1.41 GeV and 1.44 ~ 1.56 
GeV respectively, while the mass for isospin state composed of ss is 1.62 ~ 1.73 GeV. 

PACS numbers: 11.55.Hx, 14.40.Cs 
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I. INTRODUCTION 



Although the quark model has achieved great successes for several decades, the funda- 
mental structures of scalar mesons are still controversial. So far, there is not a definite 
answer on whether they are two-quark states, multiquark states or even glueball, molecule 
states among the light scalars yet . Much efforts have been given to the study of 

decay and production of these mesons. However, many theoretical predictions on properties 
of scalar mesons, in particular on the production of them in exclusive heavy flavor hadron 
decays have large uncertainties due to the complicated non-perturbative effects. 

It is no doubt that the hadronic light-cone distribution amplitudes are the important 
ingredients when applying factorization theorem to analyze these exclusive processes. The 
distribution amplitudes which are governed by the renormalization group equation can be 
obtained by integrating out the transverse momenta of quarks in hadron for hadronic wave 
functions. Unfortunately, only twist-2 light cone distribution amplitudes of scalar mesons 
have been calculated in Ref. Q in the framework of QCD sum rules jjj]. So the unknown 
twist-3 distribution amplitudes will bring obvious uncertainties to the final results. In this 
work, we investigate twist-3 distribution amplitudes of scalar mesons in order to improve 
the accuracy of theoretical predictions of the scalar mesons. 

The calculation of moments for distribution amplitudes making use of QCD sum rules 
was presented in much detail in the pioneer work of Once the moments are known, we 
can construct various models to obtain the distribution amplitudes for hadrons. Following 
the same method, we will calculate the first two non-zero moments of twist-3 distribution 
amplitudes for ao(ud), K^ds) and fo(ss) respectively based on renormalization group im- 
proved QCD sum rules. Besides, we will expand the twist-3 distribution amplitudes of scalar 
mesons according to Gegenbauer polynomials as usual and use the moments obtained to de- 
termine the first two Gegenbauer coefficients. As for and fo meson, the odd moments 
for both of the two twist-3 distribution amplitudes (see definition in Eq. (TjQ)) are zero due 
to conservation of charge parity and isospin symmetry. However, the odd moments for Kq 
meson do not vanish when including SU(3) symmetry breaking effects. 

The structure of this paper is as below: After this introduction, we derive the general sum 
rules of moments for twist-3 distribution amplitudes of scalar mesons in section [Til Then 
we will give the inputs used in our work and present the numerical results of the first two 
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moments for the above three scalar mesons in section IIHI The last section is devoted to our 
conclusions. 



II. FORMULATION 

In the valence quark model, there are two twist-3 lieht-cone distribution amplitudes for 

n 

scalar mesons which are defined as |6| 

(S(p)\q 2 (yMx)\0) = msfs t due^^^u, /*), 

Jo 

(S(p)\q 2 (y)a^ qi (x)\0) = -m s fs{p»z v -p v z,) £ due^™^^^, (1) 

with z = y — x, u = 1 — u and u being the momentum fraction carried by the q 2 quark in the 
scalar meson. Here the scalar density meson decay constant fs is defined by (S(p)\q 2 qi\0) = 
msfs- This scalar decay constant fs can be connected with the vector current decay constant 
fs which is defined as {S (p)\q 2 ^qi\0) = fsPn' 

Hsfs = fs, MS = f-^r- (2) 

with mi, m 2 and ms being the mass of q\, q 2 and scalar meson respectively. The normaliza- 
tion of these two twist-3 light cone distribution amplitudes are Jq du<j) s {u) = ducj) s {u) = 1. 
In general, they have the following form 



s (u^) = l+Y,^)C 1 J 2 (2u-l), (3) 

m=l 



jj) = 6u(l — u) 



i + E^(^ /2 (2«-i) 



(4) 



m=l 

with the Gegenbauer polynomials C\ /2 {t) = t, C 2 1/2 (t) = |(3t 2 - 1), C 4 1/2 (t) = |(35t 4 - 
30t 2 + 3), Cl'\t) = 3t, C 2 /2 (t) = f (5t 2 - 1), cf 2 (t) = f (2 It 4 - Ut 2 + 1), etc. 

Now we are ready to calculate the moments of the two twist-3 distribution amplitudes 



defined in Eq. (TTJ) making use of background field method in QCD {{], [lo|, From Eq. 
(CQ), one can easily find 

(0|ft(0)(fc. Drq 2 (0)\S(p)} =m s f s (p- z) n (C), 

ft _|_ 1 _ 

(0|fr(0)(i* D) n+1 a fJiU q 2 (0)\S (p)) = -i^m s f s { Ptl z v - Pl/ z,)(p ■ z) n (0, (5) 
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with 



du(2u- l) n <t> s s (u,fi), (O = / du(2u-ir<j>° s (u,(i). 



(6) 



In order to calculate the above scalar moments and tensor moments , we consider 
the following two different correlation functions, respectively 



• J d'xe^iOm^xXiz- g 2 (0) ?1 (0)}|0) = -(z ■ q) n I^\q 2 ), (7) 

z / d^xe^^Tiq^Wuiiz- D) n+1 q 2 (x),q 2 (0) qi (0)}\0) = i(q^ - q u z^){z ■ q) n I^°\q 2 ).(S) 



In the deep Euclidean region (— q 2 ^> 0), the correlation functions ( J71I81) can be computed 
using operator product expansion at quark level. The results with power corrections to 
operators up to dimension-six and lowest order of a s corrections are displayed as: 



7 (2n,0) (g2) 



QCD 



-ln^-(2m 1 m 2 - q 2 ) + Jil(0|G 2 |0> 



l% 2 2n + 1 /x 2 



,2n+l 



mi + m 2 )(<7i<7i) 



57T g z 
,2n+l 



1 1 

+ 2<V 



,8n + ll n . 

m 2 + n{ mi + 2m 2 

o 



m 2 + mi)(g 2 g 2 ) 



(qxaGqx) 



+ 



,8n + ll 
mi + n( m 2 + 2mi 



'~8T 
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Sn 2 - 16n - 21) 



q 2 o-Gq 2 ) | 

(gigi) 2 + fe^) 2 487ra s (gigi)(g 2 g 2 ) 
o 4 9 o 4 



(9) 



for the scalar density even moments 
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2 2\i ? i a s I 2 2^ 

2 v m x - m 2 )ln— — + ^( m i ~ m 2> 



QCD 



>7T 



1 'I , 

/i 2 



-t(0|G 2 |0) 



+ (mi + m 2 ) 



(<7i?i) - (92^2) . 107ra s (gigi) 2 - (g^) 1 
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q 

for the first moment of scalar density, and 



9s 1 f 5 5 "I 

y~ {(4 7711 + 2m 2 )(giaGgi) - (-m 2 + 2mi)(g 2 crGg 2 ) j , (10) 



QCD " 16vr 2 2n + 3 



1 + 2 In 



2 g^ 



16n 2 + 14rc + 5 



24 q 
32n 2 + 18n - 35 vra 



81 



4 [mi(gicrGgi) + m 2 (g 2 cxGg 2 )] 

1 ((<Mi> 2 + (<?2g 2 > 2 ) 



4 J 



3"7ra s 1 

"g - g 4 



(9l9l>(?2?2>5n,0, 



(11) 



for the tensor even moments, and 
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Qtt 



ml) 



2 In 



471-2 



rrir. 



1 + ln 



2 7i 



(0|G 2 |0) 
q 4 

3 miiqiaGqi) - m 2 (q2(jGq2) 
7 9s 



| m 1 (q 1 q 1 ) - m 2 (g 2 g 2 ) 27ra s (gigi) 2 - (<?2<?2) 2 



(12) 



q* 9 g 4 

for the first moment of tensor sum rule. Since we are concerned with only the first two 
Gegenbauer coefficients, we do not display the explicit forms of sum rules for other odd 
moments for simplification. When n is equal to 0, the Eq. (Q is in accord with the results 



shown in Ref. 



121 ] . On the other hand, the correlation functions fITIIHl) can also be calculated 



at hadron level by inserting a complete set of quantum states E|n)(n|, which are written as 
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lmli 2n '°\q 2 
Iml^iq 2 
lmli 2n '°\q 2 
lml^°\q 2 



had 



Mf-rni,)mimr)+* 



8vr 2 2n + 1 



(2m!m 2 -q 2 )6(q 2 -s s ), (13) 



- -nS(q 2 - r4)m%mi) + - ml)9{q 2 



'had 



-7i6{q 2 



2 -^rn 2 s f 2 s (e n )-n ' 



8vr 2 2n + 3 



q 2 6{q 2 - 



'had 



-tr% 2 - ml) 2 -mlP s (C) + ^(m 2 - m 2 )% 2 



(14) 
(15) 
(16) 



Here the quark-hadron duality has been used to obtain the above equations. Then we can 
match these two different representations of correlation functions (171181) calculated in quark 
and hadron level by the dispersion relation 

\ml is) hoo- 



ds- 



'QCD- 



(17) 



In order to suppress the contributions from the excited resonances and continuum states, 
we apply the Borel transformation to both sides of the above equation. On the other hand, 
this transformation can also remove the arbitrary polynomials in q 2 . Then we obtain 

1 1 



7TM 2 / dse ~ S/M2lmI (s)had = B M iI(q 2 ) QC D, 



(18) 



where M is the Borel parameter, and Bm 2 is the operator of Borel transformation which is 
defined as 

7 2^(n+l) / d V 



B 



m 2 



lim 

— q ,71 — >oo 

-q 2 /n=M 2 



-q 



(19) 



Finally, substituting Eqs. (l9] fT2l and (I13H16P into the Eq. (fT8l) . we have the scalar density 
even moments 



m|/|e - m |/M 2(e) 



8tt 2 2n 
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-77T.2 + 2mi 



(q 2 aGq 2 ) 
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487ra s 
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(9i9i)<9292) 



n 2 -16n-21)— ((9 ig i) 2 + (9 2 g 2 ) 



(20) 



the first moment of scalar density 
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-mlP s e-™V M *{e s ) 



8tt 2 



+ (m 2 — m 2 



(m 2 - m 2 ) / dse s/M - (mi + m 2 )((qiqi) - (9292)) 



1-ln-fW 



-(— G 2 ) + 



4M 2 7T 
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9M 2 
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,(21) 



tensor even moments 
2n + 1 



m 2 / 2 e - m |/M 2(e) 
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Jn,0 
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M 2 



(22) 



and the first moment of tensor current 



-m 2 s f 2 s e- m s/ M 2 IP 1 
3 



— -fm? — m 2 , 
4tt 2V 1 2 



/ S e~ s/M2 ds - (mi(9igi) - m 2 {q 2 q 2 )) 
Jo 



67T 



l + 21niW 



(ml -ml) M2 M2 (0\G 2 \0) 



3 mifaaGqi) - m 2 (q 2 aGq 2 ) | 27ra a faffi) 2 - (g 2 g2) 2 
V s m 2 + 9 M 2 {26) 

Here the vacuum saturation approximation [jj], Q has been used to describe the four quark 

condensate, i.e., 



= [^ADSBcSasS^SadSbc - 5 Ac^BD5 ai 5p 6 5 ac 5 bd ] (q A q A ) (q E ' q E '} . 



(24) 
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Here a,fl,j,8 are the spinor indices, a,b,c,d are the color indices, and A,B,C,D denote 
the flavor of quarks. Besides, the flavor indices in the right hand side of the above equation 
do not mean the sum of all flavors. 

It is noted that all the parameters in the above sum rules are fixed at scale of Borel mass 
M. The renormalization group equations of decay constant, quark mass and condensate are 
given as [3| 

(qq)m = (qq)?(- 

(a s G 2 ) M = (a s G 2 )^ (25) 



with b = (33 — 2n/)/3, rij is the number of active quark flavors. Making use of the orthog- 
onality of Gegenbauer polynomials 

1 dxC]l\2x - \)C]l\2x - 1) = - J— S mn , 
o 2n + 1 

/' dxx(l - x)CT{2x - l)CU\2x - 1) = (n + 2)(n + \ B , (26) 
Jo 4(2n + 3) 

the Gegenbauer moments a;, 6; can be related to moments, (^) for example: 

ai = 3(6), a 2 = ~(3<£ 2 ) - 1), a 4 = j|(35(6) - 30(6) + 3), 
h = h = ^(5(6) - 1), b 4 = ^(21(6) - 14(6) + 1). (27) 

The renormalization group equations of Gegenbauer moments are given as 

M/i)) = M^)>(^^)" 7 " /6 ' (KM) = ( 6 «(^))(^^)" 7n/6 > (28) 



where the one-loop anomalous dimensions are 



16| 



7n 5 -^(l- 7--T^r-^ + 4 E7)' 7l = qi+4E 7 , (29) 



n+l)(n + 2) ^ 2 JV 'V 



with C F = 4/3. 



III. NUMERICAL RESULTS AND DISCUSSIONS 

The input parameters used in this paper are taken as jf], Q, H, Q 
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FIG. 1: Mass (left solid line) and decay constant (right solid line) of ao from scalar sum rule in 
Eq. (|20p with ss = 4.5 GeV 2 as a function of Borel parameter M 2 . The dashed line denotes the 
contribution from continuum states in the total sum rules and the dot-dashed line is the ratio of 
dimension-six condensate contribution in the total sum rules. 
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(ss) = (0.8 ± 0.1)(uu), (uu) 9i (dd) = -(1.65 ± 0.15) x 10~ 2 GeV 3 , 

(fG a ^G a ^) = (0.005 ±0.004)GeV 4 , (g s uaGu) {gJaGd} = mg(uu), 

m u (lGeV) = 2.8MeV, (g s saGs) = (0.8 ± 0.1)(^W(rGw), (30) 

m d (lGeV) = 6.8MeV, a,(lGeV) = 0.517, 

m s (lGeV) = 142MeV, m 2 = (0.8 ± 0.2)GeV 2 . 

Here all the values for vacuum condensates are adopted at the scale fi = lGeV. Next we 
are ready to analyze the sum rules for the scalar meson nonet one by one. 

A. Mass, decay constant and moments for ao meson 

1. Determination of mass, decay constant and scalar moments (Cs^al) of ao from sum rules 
in (Op 

Here ao is the scalar meson with quark contents qiq2 = du. A common way to obtain the 
sum rules of meson mass from Eq. (1201) is taking the logarithm of both sides of this equation, 
and then applying the differential operator M A d/dM 2 to them. However, firstly we need 
to fix the value of threshold parameter and Borel parameter in order to obtain the value of 
the mass. As far as the threshold parameter ss is concerned, its value should be adopted 
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so that the Borel window is stable enough which indicates that the mass is independent 
of the choice for M 2 in some region. For the choice of the Borel parameter, one requires 
that the contribution from continuum states is less than 30% and the contribution from 
dimension-six condensates is less than 10%. In view of the above requirements, we choose 
the threshold parameter s$ = (4.5 ± 0.3) GeV 2 , such that the stable Borel window is in the 
range M 2 £ [1.60, 1.80] GeV 2 which is shown in Fig. [TJ From this figure, we can observe that 
the mass for gig2 = du scalar ground state is m ao = (1320 ~ 1410) MeV. This is very close to 
the physical state a (1450) [19|. It should be pointed out that the possibility of the existence 
of light scalar resonance near 1.4 GeV was firstly predicted by Ref. [20] as the first radial 
excitations of ao(980) according to the so called "linear dual models" on the assumptions 
of qq structure of ao(980). The decay constant of ao can be easily read from the sum rules 
in Eq. (120]) as soon as the mass is known. The decay constant within the Borel window is 
also plotted in Fig. [1] as the second diagram. It is easy to find that the decay constant is 
quite stable within the Borel window M 2 £ [1.30,1.60] GeV 2 when the contribution from 
continuum states and the dimension-six condensate is less than 30% and 10%, respectively. 
Therefore, we obtain the decay constant as / ao (lGeV) = (322 ~ 341) MeV. In the following 
subsections, all the values of decay constants and moments are calculated at scale of 1 GeV 
unless explicitly pointed out. 

From the definition of twist-3 distribution amplitudes for ao, it can be found that only even 
Gegenbauer moments are non-zero due to conservation of charge parity and isospin symmetry 
as mentioned in the introduction. Next we are going to consider the second and fourth 
moments of (p s ao from scalar density sum rules for ao meson. Just as the determination of mass 
and decay constant, one should find a stable window for the sum rule of each moment. The 
contributions of continuum states and dimension-six condensates are plotted in Fig. [2l where 
the moments (f JgJ) within the Borel window M 2 £ [1.15, 1.45] GeV 2 ([1.25, 1.55] GeV 2 ) are 
also included. For the second (fourth) moments, the contributions from both continuum 
states and the dimension-six condensates are less than 30% (35%). Then we have (£ 2 ao ) = 
0.29 ~ 0.31 and (£L ) = .16 ~ 0.19. 
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FIG. 2: (Cs,a ) (l e ft solid line) and (Cs,a ) (right solid line) from scalar sum rules in Eq. (j20j) with 
sg = 4.5 GeV 2 as a function of Borel parameter M 2 . The dashed and the dot-dashed lines are the 
ratio of contribution from continuum states and dimension-six condensates, respectively. 



2. Determination of mass, decay constant and tensor moments (£ CT ,ao ) °f a o from sum rules 
in {HP 

In the above, we have got the mass and decay constant for a meson from scalar density 
sum rules in Eq. ( BP]) . Similarly, we can also extract them from tensor sum rules in Eq. (1221) . 
Moreover, the values of mass and decay constant may not be exactly the same between these 
two sum rules due to different correlation functions adopted for them. Following the similar 
procedure, one can get the mass and decay constant: m ao = (1270 ~ 1390) MeV, f ao = 
(325 ~ 350) MeV, which are very close to the range we got from the sum rules of Eq. (j20p 
in previous subsection. The mass (decay constant) is obtained under the condition that 
the contributions from both continuum states and the dimension-six condensates should 
be less than 30% (25 %) respectively in total sum rules. The threshold parameter s s is 
still adopted as (4.5 ± 0.3) GeV 2 , while the Borel windows are M 2 e [1.60, 1.80] GeV 2 and 
[1.20, 1.50] GeV 2 , respectively. Making use of the mass and decay constant, we can determine 
the second and fourth moments (£ 2 j(I0 ), (Ctao) f° r the tensor twist-3 distribution amplitude 
of ao meson within the Borel window M 2 £ [1.20, 1.50] GeV 2 and [1.15, 1.45] GeV 2 as shown 
in Figj3l Here the contributions from continuum states and the dimension-six condensate 
are no more than 30%, which indicate that the sum rules for these two moments are reliable. 
Hence, the results for (£ 2 !ao ) and (£^ Qo ) are 0.20 ~ 0.22 and 0.093 ~ 0.12, respectively, within 
the given Borel window and threshold parameter. 
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FIG. 3: (£ 2 a() ) (left solid line) and (Cta ) (right solid line) from tensor sum rules in Eq. (|22p with 
So = 4.5 GeV 2 as a function of Borel parameter M 2 . The dashed and the dot-dashed line represent 



the ratio of contribution from the continuum states and dimension-six condensates. 



B. Mass, decay constant and moments for Kq meson 

1(2) 

1. Determination of mass, decay constant and moments (C s k*) of Kq from scalar density 
sum rules 

As explained before, here the scalar meson Kq is made up of us quarks. Different from 
the do meson, both odd and even moments of distribution amplitudes for Kq are non- 
zero. The mass and decay constant of Kq can be derived from scalar density sum rules in 
Eq. ( 1201) following the same method as done for cto case. The threshold value is chosen as 
Ss = (5.4 ± 0.3) GeV 2 in the sum rules of Eq. (1201) for Kq meson in order to gain the stable 
Borel window M 2 G [1.90,2.10] GeV 2 and [1.30,1.70] GeV 2 for mass and decay constant, 
respectively. Then we can obtain the value of mass (decay constant) of Kq as mx* = 
(1460 ~ 1560) MeV (f K * = (344 ~ 368) MeV) with the requirement that the contributions 
from both continuum states and dimension-six operator are less than 30% (25%). 

Then we try to calculate the first and second moment for Kq meson scalar twist-3 distri- 
bution amplitude according to sum rules (121]) and (1201) . respectively. For the first moment 
(£s k*) °f scalar density, we require that the contributions from both the continuum states 
and dimension-six condensates should be less than 15% in order to obtain stable Borel win- 
dow. As for the sum rules of the second moment (£ 2 K *), the contributions from both the 
continuum states and dimension-six operators are less than 20%. From the Fig. HJ we can 
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FIG. 4: {£l x*) (left solid line) from scalar sum rules (|21|) and k*) (right solid line) from sum 
rules (|20p with = 5.4 GeV 2 as a function of Borel parameter M 2 . The dashed and the dot- 
dashed lines indicate the ratio of continuum states and dimension-six condensates to the total sum 
rules, respectively. 
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FIG. 5: {CaK*) O e ft solid line) from tensor sum rules in (|23p and (S, 2 - k*) (right solid line) from 
sum rules in (|220 with Sg = 5.4 GeV 2 as a function of Borel parameter M 2 . The dashed and the 
dot-dashed lines represent the ratio corresponding to the contribution from continuum states and 
dimension-six condensates in the total sum rules, respectively. 
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read out the results of the first scalar moment K *) as (0.61 ~ 1.42) x 10 2 within the Borel 
window M 2 E [1.90,2.20] GeV 2 , and the second moment as 0.29 ~ 0.33 within the 

Borel window [1.20, 1.50] GeV 2 . The threshold parameter is fixed at ss = (5.4 ±0.3) GeV 2 . 
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2. Determination of mass, decay constant and moments (C a ) of Kq from tensor current 
sum rules 

Similarly, we can also derive the results of mass and decay constant from the ten- 
sor operator sum Rules in Eq. fl22|) . Here we will only show our values of them as 
m K , = (1440 ~ 1550) MeV and f K * = (349 ~ 375) MeV within the Borel window 
M 2 e [2.00,2.20] GeV 2 and [1.30,1.60] GeV 2 . The threshold parameter is set the same 
as before, s s = (5.4 ± 0.3) GeV 2 . The contributions from both the continuum states and 
dimension-six condensates are required to be less than 30% (20%) for mass (decay constant) 
sum rules respectively. The mass of Kq determined here are consistent with the one deter- 
mined in the previous subsection from sum rules (|20|) . which is quite close to the physical 
state i^o(1430). The first and second moments (£,a,K*) of the tensor twist-3 distri- 

bution amplitude can be computed following the same method, which have been plotted in 
Figj5j It is found that the sum rules for these two moments are quite stable within the Borel 
window M 2 e [2.20,2.60] GeV 2 and [1.00,1.20] GeV 2 , respectively, since contributions from 
both continuum states and dimension-six condensates are less than 10%. The results for 
them are shown as: = (2.2 ~ 3.3) x 10~ 2 , (& >K *) = 0.20 ~ 0.25. 

C. Mass, decay constant and moments for fo meson 

1. Determination of mass, decay constants and scalar moments of fo from Sum Rules 

in (MW 

Here fo refers to the scalar meson which is made up of ss quark. The sum rules for f 
are much the same as for the ao meson. The odd moments vanish due to conservation of 
C parity. Therefore, we will only consider the first two even moments, (£ 2 / ) and (^/ ) 
of the scalar twist-3 distribution amplitude for the fo meson. Since the calculations are 
similar as that for ao meson, the mass and decay constant for fo are given straightforward 
as m fo = (1640 ~ 1730) MeV and f fo = (369 ~ 391) MeV within the Borel window 
M 2 G [2.50,2.70] GeV 2 , [1.70,2.00] GeV 2 respectively. The threshold parameter is set as 
Ss = (6.5±0.3) GeV 2 . Here we also require that the contributions from both the continuum 
states and dimension-six condensates are less than 30% (20%) for mass (decay constant) 
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FIG. 6: (£s / ) (left solid line) and (Ctj ) (right solid line) from scalar sum rules in Eq. ({20]) with 
ss = 6.5 GeV 2 as a function of Borel parameter M 2 . The dashed and the dot-dashed lines reflect 
the ratio of continuum states and dimension-six condensates to the total sum rules, respectively. 

sum rules respectively. As for the second and forth moments (£ 2 / () ),(£s t ), the results within 
the same Borel window M 2 e [1.60, 1.90] GeV 2 are plotted in Fig. [61 The number of (£ 2 / ), 
(£a f ) are 0.29 ~ 0.31 and 0.17 ~ 0.20 respectively within the given Borel window and 
threshold parameter. The requirement that the contributions from the continuum states 
and dimension-six operator are less than 25% (30%) for the second (forth) scalar moments 
has been used. 

2. Determination of mass, decay constant and tensor moments (C^fl) °f fo from Sum Rules 
in (dj> 

The mass and decay constants of fo can also be derived from tensor sum Rules in Eq. (1221) . 
Adopting the same threshold parameter as the scalar density sum rules, we obtain the 
results as m fo = (1620 ~ 1710)MeV and f fo = (381 ~ 426)MeV within the Borel window 
M 2 G [2.50,2.70] GeV 2 , [1.20,1.60] GeV 2 respectively. Here we require that contributions 
from both continuum states and dimension-six operators are less than 30% (10%) for mass 
(decay constant) sum rules. The mass we get here from tensor sum rules and also that from 
the scalar density sum rules (127)1) in previous subsection is close to the physical state /o(1710). 
The second and forth moment (£ 2 f Q ), (^/ ) °f tensor twist-3 distribution amplitude are also 
displayed in Fig. within the Borel window M 2 e [1.50,1.80] GeV 2 and [1.60,1.90] GeV 2 
respectively. The condition that the contributions from the continuum states and dimension- 
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FIG. 7: (Cct/ ) (l e ft solid line) and f ) (right solid line) from tensor sum rules in Eq. (|22|) with 
sZ = 6.5 GeV 2 as a function of Borel parameter M 2 . The dashed and the dot-dashed line represent 



the ratio of contribution from the continuum states and dimension-six condensates. 



TABLE I: Masses, decay constants and Gegenbauer moments from the scalar density sum rules 
(I20|2ip at the scale fi = lGeV and 2.1GeV (shown in the second line of each meson) 



state 


m(MeV) 


/(MeV) 


ai(xl0~ 2 ) 


«2 


04 


a 


1320 ~ 1410 


322 ~ 341 
391 ~ 414 





-0.33 0.18 

-0.26 ~ -0.14 


-0.11 ~ 0.39 
-0.075 ~ 0.27 


K* 


1460 ~ 1560 


344 ~ 368 
418 ~ 447 


1.8 ~ 4.2 
1.6 ~ 3.8 


-0.33 - -0.025 
-0.26 ~ -0.020 




fo 


1640 ~ 1730 


369 ~ 391 
448 ~ 475 





-0.33 ~ -0.18 
-0.26 ~ -0.14 


0.28 - 0.79 
0.19 ~ 0.54 



six operators are less than 25% (30%) is adopted for the second (forth) tensor moment. The 
value of (£ct/ ) an< ^ (£tf ) are ~ an d 0.057 ~ 0.082 within the given Borel window 
and threshold parameter. 

Now we have finished the calculation of the moments (O^) of twist-3 distribution ampli- 
tudes for scalar mesons ao, Kq and /o in the framework of QCD sum rules. With the results 
of (Q a ))j ^ i s straightforward to derive the Gegenbauer moments a m and b m in Eq. ( 13141) 
using Eq. (!27j) . The results for the first non-zero Gegenbauer moments at lGeV and 2.1GeV 
scales are shown in table [J and [TTJ They can be applied to various approaches involving light 
cone distribution amplitudes of hadrons, such as perturbative QCD approach 21| . QCD 
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TABLE II: Masses, decay constants and Gegenbauer moments from the tensor sum rules (|22I23|) 
at the scale /i = lGeV and 2.1 GeV (shown in the second line of each meson) 



state 


m(MeV) 


/(MeV) 


6i(xl0- 2 ) 


b 2 


6 4 


a 


1270 ~ 1390 


325 ~ 350 
395 ~ 425 





~ 0.058 
~ 0.041 


0.070 ~ 0.20 
0.045 ~ 0.13 


K* 


1440 ~ 1550 


349 ~ 375 
424 ~ 456 


3.7 ~ 5.5 

2.8 ~ 4.2 


~ 0.15 
~ 0.11 




fo 


1620 ~ 1710 


381 ~ 426 
463 ~ 518 





-0.15 ~ -0.088 
-0.11 0.062 


0.044 ~ 0.16 
0.028 ~ 0.10 



factorization approach 



231 ] etc. As mentioned above, the odd 



and light-cone sum rules 
moments of twist-3 distribution amplitudes for scalar mesons and /o are zero due to 
conservation of charge parity and flavor symmetry as explained in the introduction. As a 
byproduct, we also collect the masses and decay constants of scalar mesons in table HI and HI 
These masses indicate that the ground state of qq scalars are probably a (1450), 7^(1430) 
and / (171Q1. 

In ref. 24J, the authors also studied the mass and decay constant of scalar meson Kq. 
Their results are m^* = (1410 it 49) MeV and fx* = (427 ± 85) MeV, which are consistent 
with our results within error bar. 



IV. SUMMARY 



In this work, we have studied the masses, decay constants and twist-3 distribution am- 
plitudes of scalar mesons based on the renormalization group improved QCD sum rules. It 
is shown that the mass sum rules for scalar mesons are not very satisfied, since the Borel 
windows are a bit narrow for all the three scalar mesons. Our results for the scalar meson 
masses show that the physical states a (1450), _K"q(1430) and / (1710) are preferred to be 
the ground state of scalar mesons. The sum rules for decay constants of these three scalar 
mesons are very stable in a much broader Borel window. The second and forth scalar mo- 
ments of ao can be obtained with 30% and 35% uncertainties respectively, while both the 



1(3 



second and forth tensor moments of ao can be derived within 30% uncertainties. As for the 
Kq meson case, the first and second moments of scalar density twist-3 distribution ampli- 
tude 4> S K * are obtained under 15% and 20% uncertainties respectively. The uncertainties can 
be reduced to 10% for both the results of the first and second moments of tensor twist-3 
distribution amplitude 0^* . For the case of /o meson, the second moment for both of scalar 
twist-3 distribution amplitude 0j o and tensor twist-3 distribution amplitude 0j o each has 
25% uncertainties. Besides, the fourth moment for each of these two distribution amplitudes 
could be obtained within 30% uncertainties. It is also worthwhile to emphasize that the cor- 
relation functions are calculated to leading a s power based on operator product expansion 
in this work, which will bring some additional uncertainties to mass, decay constants and 
Gegenbauer coefficients. 

It is found that the second Gegenbauer coefficients of scalar density twist-3 distribution 
amplitudes for Kq and fo mesons are quite close to that for ao, which indicates that the 
SU(3) symmetry breaking effect is tiny here. However, this effect could not be neglected for 
the forth Gegenbauer coefficients of scalar twist-3 distribution amplitudes between ao and fo- 
Furthermore, one can also observe that the first two Gegenbauer coefficients corresponding 
to tensor current twist-3 distribution amplitudes for all the ao, Kq and fo are very small. 
As is well known, the light-cone distribution amplitudes play a critical role for hadronic 
decay processes in the framework of factorization theorem where it describes the bound 
state effect of hadrons. The available twist-3 distribution amplitudes of scalar mesons allow 
us to improve the accuracy of the theoretical predictions on the properties of scalar mesons, 
in particular for the heavy flavor hadron decays to scalar mesons; so that it is very helpful 
for us to understand the structure of scalar mesons and strong interactions. 
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